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The gravitational radiation degrees of freedom of freedom 
are described in the framework of the 3+1 decomposition of 
spacetime. The relationship with eigenfields of the Kidder- 
Scheel-Teukolsky (KST) equations jjj is established. This re- 
lationship is used to fix a parameter in the KST equations 
which is related to the ordering ambiguity of space derivatives 
in the Ricci tensor, which is inherent to first order evolution 
systems, like the ones currently used in Numerical Relativity 
applications. 



I. INTRODUCTION 

The structure of Einstein field equations has deserved 
great interest since the very beginning of General Rel- 
ativity. It was early noticed that, by rearranging the 
order or partial derivatives, the principal part of the four 
dimensional Ricci tensor could be written as a sort of 
generalized wave equation 0,0] 



2R IU , = -Og^ + B^ + duT,, 



(1) 



where the box stands for the d'Alembert operator acting 
on functions and we have written for short 



(2) 



This opened the way to the use of harmonic spacetime 
coordinates (Dx^ = 0) in order to obtain an hyperbolic 
evolution system |^-^| . 

By the middle of the past century, however, the inter- 
est focused in the relativistic Cauchy problem. The 3+1 
decomposition of the line clement: 

ds 2 = -a 2 dt 2 

+ 7^ (dx i + ft dt) (dx j + f3 j dt) i,j = 1, 2, 3 (3) 

allowed one to express six of the ten second order original 
equations as a system of evolution equations for the met- 



ric 7^ 



and the extrinsic curvature Ka of the t = constant 



slices, namely: 

(d t - £fs)lij = -2 a Kij 
(d t - £>p)Kij = —ViOij + a 



(4) 



2Kf 



trKKij] (5) 



where C stands for the Lie derivative (we restrict our- 
selves to the vacuum case for simplicity) . The remaining 
four equations could instead be expressed as constraints: 



{3) R-tr(K 2 ) + (trK) 2 = 
V fe K k i-di(trK) = 



(0) 
(7) 



This opened the door to a new way of obtaining hyper- 
bolic evolution systems The key point was to use 
in one or another way the momentum constraint to 
ensure hyperbolicity while keeping the freedom of choos- 
ing arbitrary space coordinates on every t = constant 
slice. This allows for instance to use normal coordinates 
(f3 l = 0) without affecting the mathematical structure of 
the evolution system. (See Ref. jHj] for a detailed compar- 
ison of the "old" and "new" hyperbolic systems). 

These findings came at the right moment for people 
working on Numerical Relativity with a view on the grav- 
itational waves detector projects starting by the turn of 
the century. Following the wake of these first works, 
many groups found their own way of combining the mo- 
mentum constraint with the evolution equation (^J), lead- 
ing in each case to a new brand of hyperbolic systems 



Suddenly, the problem of having too few hyperbolic 
formalisms for Numerical Relativity turned into the op- 
posite problem of having too many of them. Some of the 
works considered even multiparameter families of hyper- 
bolic systems |I| Jlo| , [L4| . Faced with the problem of choice, 
we think that the right question at this point is: what 
are hyperbolic systems for? 

There are many answers, of course, but we can get a 
hint when we realise that hyperbolic systems basically de- 
scribe propagation phenomena along characteristic lines. 
This sets up the question of whether we can use these 
formalisms in the context of gravitational wave propa- 
gation. The 3+1 formalism does not seem at first sight 
very useful in that context, because it describes space- 
time as a foliation of t = constant hypersurfaces. Then, 
the natural structure associated with gravitational waves 
in that formalism is not the light cone itself, but the two- 
dimensional surfaces (wavefronts) obtained as the inter- 
section of the light cone with the t = constant hyper- 
surfaces. The geometry of wavefronts will be discussed 
in Section II. The second difficulty arises from the fact 
that, because of the constraint equations (^, [?]), only 
two degrees of freedom correspond to gravitational ra- 
diation in spite of the complete (strong) hyperbolicity 
of the formalisms. This point can be easily illustrated 
by looking at the list of characteristic speeds given in 
the KST paper M namely {0, ±1, ±ci, ±C2, ±03} where 
ci, C2 and C3 depend of a list of arbitrary parameters 
related with the gauge choice (ci) and with the way in 
which the constraints are used to enforce hyperbolicity. 
It is clear that the degrees of freedom propagating with 
speed 0, ±ci, ±C2 or ±03 can not describe gravitational 
radiation, which we know that propagates at light speed 
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(v = ±1), independently of how we manage to write down 
our equations. The two remaining degrees of freedom, 
as expected, are related with transverse traceless com- 
ponents of the extrinsic curvature Kij. We will analyse 
more closely these eigenfields in Section III. 



II. LIGHT CONES AND WAVEFRONTS 

Let us consider the propagation of a burst of gravita- 
tional waves. From the geometrical point of view, it can 
be described as a foliation of spacetime by a set of null 
hypersurfaces (light cones), namely 

<fr(x, t) = constant (8) 

d4-d4 = Q (9) 

From every t — constant hypersurface, the propagation 
is seen instead as a succession of wavefronts. It can be 
obtained by setting t = t a in (||) to obtain 

4>(x, t ) = constant (10) 

so that the null foliation (||) induces a spacelike foliation 
( |l0| ) by two-dimensional surfaces (wavefronts) on every 
t = constant slice. 

Let us consider now an adapted coordinate system on 
the t — t Q slice so that the space unit normal n k to the 
wavefront (|l0| ) is given by 

n k = N 5 k z (11) 

and the other two coordinates x a (a — 1,2) display the 
wavefronts surface. The line element for the t — t Q slice 
can then be written as (2+1 decomposition) 

7y dx i dx j = N 2 dz 1 

+ a ab {dx a + \ a dz) (dx b + \ b dz) a, b= 1,2 (12) 

where a a b is the induced metric on every wavefront sur- 
face, so that it describes the intrinsic geometry of the 
wavefronts. 

The extrinsic geometry of the wavefront foliation jl0| ) 
is in turn described by the second fundamental form n a b, 
namely 



K a b = {dz ~ £\)o- a b 



(13) 



Both <j a b and n ao are 2+1 tensors, in the sense that they 
transform in a covariant way under changes of coordi- 
nates of the form 

z = f(z) (14) 
Z a = h a (x,z) (15) 

which transform the wavefront foliation (|To|) into itself. 

Notice that o~ a b, K a b can be easily expressed in terms of 
the three-dimensional metric 7^ and the spacelike unit 
normal n k , namely 

Vab = lab (16) 
Kab = ~n k T k ab (17) 

where T k ij stand for the Christoffel symbols of 7^. 



III. SOLVING THE ORDERING AMBIGUITY 

Let us come back now to the results of the KST paper 
|l| . The principal part of the first order evolution systems 
presented there can be written in Flux-Conservative form 
for the array of variables u = {7^, Kij, d k ij}, namely 



-d t u + d k F k {u) 



(18) 



so that propagation along the direction given by n k can 
be studied by solving the following characteristic eigen- 
value problem 

F n {u)=n k F k {u) = vu (19) 
From the transverse components K a b, one has |l| 



F n (K ab ) = ~[d n ab -(l + {) d (ab) T 



m -f a b (d nk k 



d k kn 



)] 



(20) 



so that the parameter m (corresponding to 7 in Q]) does 
not contribute to the transverse traceless part. On the 
other hand, one has [Qj 



F n (d nab ) 

F n (d a bn) 



2 K ab 

_V 
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X lab (K nn - tvK) (21) 
lab (K nn - trK) (22) 



so that it is clear again that the parameters ry, \ do not 
contribute to the transverse traceless part. 

This means that the traceless part of the combinations 



K a b ± \ Kab - (1 + C) d (ab) n ] 



(23) 



corresponds to eigenfields propagating with light speed 
(v = ±1, respectively). The parameter £ appearing in 
( p3| ) is related with the ordering ambiguity inherent to 
first order formalisms, where the space derivatives d k ij 
are considered to be independent quantities, so that the 
identity 



d[rd s ]ij = 



(24) 



is not granted and must be considered as a supplemen- 
tary constraint 0]. The value £ = — 1 corresponds to 
a decomposition of the three-dimensional Ricci tensor of 
the form (|l|), which is closer to the wave equation. The 
value C = +1 corresponds instead to the simpler decom- 
position 



1 j — d k r k ij — u(ji 



d ( ,T k n 



r kr r ij r r .i r k j (25) 



Notice that, allowing for (|17|), the combinations (]23|) in 
the C = +1 case can be written as 



K a b ± Kab 



(26) 



so that they have a geometrical meaning in terms of the 
extrinsic curvature of the t = constant hypersurfaces and 
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the wavefront foliation. For any other value of £, how- 
ever, the eigenfields described by the traceless part of (|2^ ) 
are not covariant under the coordinate transformations 
( |l4| ) and do not admit then any consistent, coordinate 
independent, physical interpretation. It follows that one 
can solve the ordering ambiguity for the derivatives in the 
principal part of the Ricci tensor by requiring that the 
eigenfields corresponding to the transverse traceless de- 
grees of freedom do have a clear geometrical meaning in 
terms of the geometry of the t — constant hypersurfaces 
and that of the wavefronts. This requirement amounts 
to select the particular value £ = +1, corresponding to 
the classical decomposition (pa). 



IV. CONCLUDING REMARKS 

The parameter dependence of the eigenfields ( p3| ) may 
seem surprising at the first sight, because these trans- 
verse traceless degrees of freedom are related with a real 
physical phenomenon (gravitational waves), which is in- 
dependent on how we do actually decompose the Ricci 
tensor. This paradox can be solved by noticing that the 
Fluxes appear in the evolution equations ([l8|) under a di- 
vergence operator; the parameter C appears when includ- 
ing the "rotational" constraint (B4j), but all these terms 
just cancel under the divergence operator, so that all val- 
ues of £ are physically equivalent. 

The special choice £ = +1 corresponds to the only 
case in which the Fluxes themselves inherit the covari- 
ance properties of the equations. We claim that this is 
an important feature, because then gravitational radia- 
tion features are explicitly inherited by the eigenfields of 
the system. 

On the other hand, there are indications that £ = +1 
may not be convenient from the Numerical Relativity 
point of view. Most of the formalisms that are being cur- 
rently (and successfully) used to that end actually imply 
the opposite choice (£ = —1). These include at least the 
Bona-Masso ||Jl^], Shibata-Nakamura pq ], Baumgarte- 
Shapiro [^7j and Anderson- York formalisms. The 
reason may be that with the ( = — 1 choice the transverse 
traceless part of the evolution system looks like that of 
the wave equation, as pointed out in jlj , which is known 
to admit a symmetric-hyperbolic form. 

Symmetric hyperbolicity is a stronger requirement 
than just strong hyperbolicity in order to ensure exis- 
tence, unicity and stability of the solutions arising from 
a given set of initial data (well posedness) When 
the evolution system is just strongly hyperbolic, but not 
symmetric hyperbolic, then well posedness is not ensured 
for non-smooth data, like the ones arising in numerical 
experiments due to finite machine accuracy. 

We actually have some preliminary evidence that the 
choice £ = +1 is prone to numerical instabilities. To con- 
firm this, one must proceed to a systematic exploration of 
parameter space like the one presented in the KST paper 



|l[ for the ( = — 1 case. There is also some theoretical ev- 
idence that the choice £ = +1, when combined with 
conventional choices of the other parameters introduced 
in Q, docs not lead to a symmetric hyperbolic system. 
This opens the way, from the theoretical point of view, 
to using the symmetric hyperbolicity requirement, in the 
C = +1 case, to restrict the values of the other parame- 
ters, then reducing the volume of parameter space. 

We are currently working along these lines, both from 
the theoretical and from the numerical point of view. 
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